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Abstract 

In this paper we present a feasibility study for measuring the combination of 
the Unitarity Triangle angle 2/3 + 7 with a time dependent Dalitz analysis in B° 
decays using D~K°ir + final states following the method proposed in pp. For this 
study we build a model for this decay using the latest experimental information and 
we investigate the possibility of fitting together with 2/3 + 7 an the relevant strong 
amplitudes and phases of the decay model. 

1 Introduction 

The interference effects between the b —>■ c and b — > u decay amplitudes in the time- 
dependent asymmetries of B decaying into D^it and D^p final states allow for the 
determination of 2/3+7. Nevertheless the extraction of 2/3+7 depends upon the knowl- 
edge of the absolute value of the ratio (r) of the amplitudes ( e.g.: r = \A(B — > 
D~7T + )\/\A(B° — > D~7i + )\) and of the strong phase difference, which nave not been de- 
termined experimentally so far. Furthermore, the ratio r is expected to be rather small, 
of the order of A| V^/V^l ~ 0.02. In addition the extraction of 2/3 + 7 from the two-body 
decays suffers from a eight-fold ambiguity. One can estimate the ratio r by deducing for 
instance the branching fraction Br(B° — > D^~ir + ) from the measured branching frac- 
tion Br(B° — > Dg*^ 7r + ). This approach is valid up to a theoretical uncertainty related to 
SU(3) breaking effects and to the size of annihilation processes which do not contribute to 
Dg tt + final states. Using all the experimental constraints 0,0, a tentative extraction 
of 2/3 + 7 nas been made, attributing 100% uncertainty to these assumptions, and getting 
[S] : 2/3 + 7 — (±90 ±46)° (the ± in front of the central value refers to the fact that 2/3 + 7 
is determined up to a 7r ambiguity). 

Another way of measuring 2/3 + 7 is to perform a time dependent Dalitz analysis 
of the three-body decay : B° — > D~ K°tt +1 , as proposed in pQ. This method allows 
in principle to overcome the disadvantages encountered in the determination of 2/3 + 7 

1 Charge conjugation is implied throughout the paper unless otherwise stated 



from D^tt or D^*'p decays. First of all, in the region of the Dalitz plane where most 
of the interference takes place the value of the ratio r is expected to be of the order of 
0.4 since it involves b — > c and b — > u colour suppressed diagrams. Furthermore 2(3 + 7 
can be measured with only a two- fold ambiguity as explained in In addition, this 
method could be theoretically clean, because the strong amplitudes and phases can be, 
in principle, measured from data. 

Early studies can be found in jlj. In this paper we present a feasibility study of the 
analysis under realistic conditions. Unlike to previous studies, we have used a model 
for the Dalitz structure of the B° — > D~K°7r + decays that uses the latest experimental 
measurements. We have also included realistic B° — B° tagging performances and explored 
the possibility of fitting the strong amplitudes and phases of the most relevant interfering 
resonances by using both the tagged and untagged events. Finally we have studied the 
effect of background in the determination of 2(3 + 7. 

2 The B° L>-i^°7r + decays 

2.1 Time dependence : general case 

The measurement of 2(3 + 7 can be achieved through a study of the time-dependent 
evolution of B° decay sensitive to that phase. This can only be performed on tagged 
events, that means events for which the original flavour of the reconstructed B is known. 
In a B factory this is possible through the determination of the flavour of the other B in 
the event since the two B mesons are produced in a coherent J PC = 1 state. Thus if 
we determine that at a time t we have a B°, then the other B in the event at the same 
time t is a B°. In practice, one of the B (the -Stag) is reconstructed in a flavour tagging 
state. The time-dependent decay rates taking into account CP violation for a B° (P+) or 
B° (P-) tagging meson are given by : 

P± = — e~ t/T (l =F C cos(Am d t) ± S sin(Am^)) (1) 
4r 

where r is the B° lifetime, Am d is the B° —B° mixing frequency and N a normalisation 
factor. The parameter S is non-zero if there is mixing-induced CP violation, while a 
non-zero value for C would indicate direct CP violation. 

2.2 The B° -> L> if°7r+case 

The final state D~K°ir + can be reached through the diagrams shown in Figure ^ and 
considering the B° — B° mixing. The parameters S and C of Eq. ^ will depend on the 
position in the B Dalitz plot. Indeed, a given bin in the Dalitz plot gets contributions 
from amplitudes and strong phases of b — > c and b — > u transitions. 

The model assumed for the decay parametrises the amplitude A at each point k of 
the Dalitz plot as a sum of two-body decay matrix elements and a non-resonant term 
according to the following expression : 

4*(« fc )e w «*<«*> = Yl a^BWlim, T, s) + a nr e 1 ^ (2) 
3 




Figure 1: Feynman diagrams describing the processes contributing to the decays B° — > 
D-K°tt+. a) B° -»■ D-K* + (K°n+) and higher K** resonances, b) B° -> D**°(D-tt + )K° , 
c) B° -> D s ** + (D + K°)ti-, d) B° -> D**°(D + 7i-)K°. The processes in a,b (c,d) are V cb 
(V u b) mediated. 

where indicates the Cabibbo allowed (suppressed) decay in each point k of the 

Dalitz plot. Each term of the sum is parametrised with an amplitude (aj or a nr ) and 
a phase (Sj or 4> nr ). The factor BWl(m,T, s) gives the Lorentz invariant expression for 
the matrix element of a resonance j as a function of the position k in the B Dalitz plot; 
the functional dependence varies with the spin s of the resonance according to the isobar 
model [H] . The total phase and amplitude are arbitrary. We have chosen amplitude unity 
and phase zero for the mode K* + (892) decaying into K°tt + . 

Four transitions have to be considered associated to probabilities of connecting a B° or 
B° initial state to a D + or D~ final states. They are indicated in Table 
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Table 1: Summary of the amplitudes involved in the 2(3 + 7 analysis. The index k refers 
to the position in the Dalitz plot. The third column is derived from the second one using 
CP transformation. 

The time dependent evolution of these probabilities can be obtained from the resolu- 
tion of the Schroedinger equation. This leads to a new formulation of Eq. ^ in which the 
first argument within the parentheses refers to the production state of the B while the 
second one indicates the reconstructed final state 2 . 



2 If one neglects b — > u transitions, the charges of the particles in the final state tag the B decay 
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P(B°, D + K\-) = Ck u « e - rBt {l - C k cos(Am^) + S 1 ; sm(Am d t)} (3) 

_ A 2 + A 2 

P{B\D + K\-) = Ck - Uk e- rBt {l + C k cos{Am d t) - S{sin(A m(i t)} 
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(4) 



with : 



fc _ A l k ~ A l k . ck _ 2Im(A Ck A Uk e^^- s ^) , cfe _ 2Im(A Ck A Uk e^^-^) 



(5) 

Using these relations, and because of the presence of the terms BWl(m,T, s) which 
vary over the Dalitz plot, we can fit the amplitudes (aj) and the phases (8j) of Eq. |2l 
together with 2(3 + 7 with only a two-fold ambiguity. 



2.3 The 5 U D~K\ + model 

We now discuss the hadronic model for the B° — > Z^.ff'V+decay. The resonances taken 
into account are listed in Table El for both V c b and V u b mediated modes and shown in 
Figure 121 The Dalitz plot can be modelled in terms of the following resonances : 

D~ — X X -f K°tt + (fsT*(892) ± , J RT*(1430) ± , i^*(1430) ± , iT*(1680) ± ) (6) 

K°-Y Y -> L)-7T+ (L)*(2400)°,^(2460)°) 
7T - Z Z -> D-JT (D Sj2 (2573) ± ) 

The i^°7r + resonances (If* like) can only come from V c b mediated processes, while both 
V c b and V u b transitions contribute to D~n + resonances (D** like). Finally the D~K° 
resonances (_D**-like) come only through V u b mediated processes. For the final state 
under study, "non resonant" contributions could come from higher K excited states or 
from higher excited D** states. In the first case only V c b processes can contribute while 
in the second both V c b and V u b processes are contributing. 

The strong phases are not known experimentally and have to be chosen arbitrarily. 
For the values of the amplitudes we can use some information available on branching 
fractions |Zj,(2]- So far we have some information which can help defining the model for 
the b — > c part. The Dalitz plot of this decay has been partially measured in jH], 0. Two 
measurements are available to partially define the D** part of the Dalitz plot : 

Br(B° -> D-R* + ) = (4.6 ± 0.6 ± 0.5) x 10~ 4 (7) 
Br(B° -> D**- 2 K + ) x Br{D*j% -> D°n-) = (1.8 ± 0.4 ± 0.3) x 10~ 5 

flavour and one is left with the standard mixing formulae : P(B a ,B°) = P(B°,B°) oc (1 — cos(Arridt)) 
and P(B°, B°) = P(B°, B°) oc (1 + cos(Am d t)) 




Figure 2: Various resonances contributing to the B° —> D K°n + Dalitz plot. From top to 
bottom, left to right : fsT*(892) ± , iT*(1430) ± , fsT*(1430) ± , fsT*(1680) ± , £>*(2400) , D*(2460)° 
and D Si2 (2573) ± 



The production of the charged D** through the process indicated in the second equation 
(Eq. |HJ) involves only V c h transitions with diagrams at the tree level (T) (similar to the 
diagram in Figure ^ (a), with a D*f~ 2 produced in the lower part and a K + emitted from 
the W), while we are interested in the neutral D** production which proceed via colour 
suppressed decays (C) (see Figure HJ-(b)). For our numerical evaluation of the branching 
fractions of interest for this analysis we have considered that the ratio \C/T\ could vary 
between 0.3 and 0.5. Nevertheless we are conscious that it is not straightforward to 
link production of specific hadronic states mediated by T and C suppressed processes. 
Furthermore, only 2 + final states have been measured so far. We make the hypothesis 
that the + states decaying into Dir modes are as abundant as the 2 + states. The possible 
existence of a non resonant contribution is considered below. 

For the b — > u counterpart we use the hypothesis that r is equal to 0.4 (where r in 
this case is defined as r = ^oT^lo^o )- The value of 0.4 is motivated by the fact that 

r = ofa x = x The term + = 0.408 ± 0.016 

0. The term ^„ ^ | can be expressed as the ratio of two colour suppressed diagrams 
involving and V uh processes. This ratio is expected to be around unity with an error 
which is difficult to estimate. For studying the impact of the r parameter on 2(3 + 7 
sensitivity we consider a variation for this parameter in a relatively wide range : between 
0.3-0.5. 

We also use the measurement of the following ratio 
B° ~D-*-+(892)(*V) 



B° -> D~K° 



7T 



to distribute the rest of events in the Dalitz plot. It has to be considered that 
no information on how to distribute these events between different excited K states 
(K*(U30), #2(1430), #*(1680)) is available. The contribution from the D** resonance 
in b — > u transitions is difficult to evaluate. By analogy we can consider a similar decay 
mode : B° — > D* + D*~ K° . These decays are mediated by tree diagrams, however the con- 
tributions from B° — > D^D^ final states are small [TO] - This result would indicate the 
predominance of non resonant or very large states. A conservative choice has been made 
not to include them in the present model. The contributions from the different resonances 
are fixed at the values given in Table 121 Finally we make the conservative hypothesis to 
have, all over the Dalitz, a small non resonant contribution and in agreement with the 
Standard Model we assume 2/3 + 7 = 2 radians (see for example [5]). All these parameters 
will be varied to evaluate the impact of the chosen model on the sensitivity of 2/3 + 7. 
The number of -B — > D~K°n + signal events is estimated using the measured branching 
fraction and observed yield jH]. We have generated 250 signal events per unit of 100 fb~ . 



2.4 Sensitivity study for B » — » D K 7r + decays 

In order to show which are the regions in the Dalitz plot that contribute most in the 
determination of 2/3 + 7, we have evaluated, on an event by event basis, the second 
derivative with respect to 2/3 + 7 of the log-likelihood ^(2^+7 ) constructed according 
to Eq0] (this likelihood is signal only and does not consider resolution effects). A very 
high statistics Monte-Carlo sample of signal events has been generated according to the 
nominal model described in Section 12.31 The result is shown in Figure El where each 
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Table 2: Lzst of mass, widths and quantum numbers of the resonances considered in our 
model, as taken from PDG200A. The last four columns present the chosen values of the 
coefficients aj and Sj in Eq. [H for the Cabibbo allowed and Cabibbo suppressed decays 
respectively. Note that the choices for the phases are arbitrary and are not indicated in 
the Table. In the numerical analysis we have evaluated the effect of a different choice of 
strong phases. 



event is weighted by the value of the second derivative with respect to 2/5 + 7 of the 
log-likelihood: 

. , d 2 lnL . , 

weighi = ¥WT7) (9) 

The choice of this weight is motivated by the fact that the uncertainty on 2(3 + 7 is 

equal to / g ilnL . Figure El shows the distribution of the Monte Carlo events in 

y ^ e 2 (2/3+7) 

the Mp +K0 versus M^ 07r _ plane (black dots) and the corresponding sensitivity weights 
superimposed (coloured squares). The regions with interference between B° — ► D**°K° 
and B° — > D**°K° colour suppressed processes show the greater sensitivity to 2/5 + 7. A 
particularly sensitive zone is at the intersection between B° — ► D~K* + and the colour 
suppressed B° — > D**°K°. In the same figure we can also notice that some sensitivity is 
also present because of the interference of the _D Sj 2(2573) ± with the excited K resonances 
(A" *(1430) ± , iT 2 *(1430) ± , iT*(1680) ± ). As expected, we observe no sensitivity in the regions 
where there is only one path (b — > c or b — ► u) to reach the DK°ir final state and a maximal 
sensitivity when there is an overlap between resonances due to both b — ► c and b — ► u 
transitions. 



3 Feasibility study 

3.1 Dependence on the Dalitz model of the 2/3 + 7 error. 

The first study is to consider the dependence on the determination of 2/3 + 7 from the 
Dalitz structure of the decay model. For this study we have used the time evolution 
equations defined in Eq. 0] without including the effect of time resolution. Preliminary 
studies using a realistic resolution indicate an increase on 2/5 + 7 uncertainty of about 
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Figure 3: a) Dalitz distribution of the very high statistics Monte-Carlo sample of signal 
events. Each event is entering into the plot with a weight given by the value of the second 
derivative with respect to 2(3 + 7 of the log-likelihood. The black points correspond to the 
same events with weight equal to unity, b) Distribution similar to Figure-a except that we 
have fixed a maximum value for the weight to be plotted in order to see in a finer way the 
structure of the weights over the Dalitz. 
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Figure 4: Average absolute uncertainty obtained on the parameter 2/5 + 7 as a function 
of the integrated luminosity for different models. The magenta thick stars refers to the 
model in tabled while the thin stars are obtained by varying the strong phases values; the 
red (dark) full circles and the green (clear) open square are obtained multiplying the D**° 
amplitudes by a factor 1.3 and 0.7 respectively. The red full (dark) and green (clear) open 
triangles refer to previous models considering r=0.5 and r=0.3 respectively. Finally the 
blue crosses refers to the model in table\^where the non-resonant component is set to zero. 
The cyan (clear) squared are just given as reference and correspond to the simplified model 
presented in (4( in which all V u b processes are considered through non-resonant amplitudes. 

10%. Instead, we have included realistic B° — B° tagging performances in terms of purities 
and efficiencies. The effect of the background is studied in the following section. For this 
study we just fit 2/3 + 7 fixing all the other parameters describing the Dalitz decay model. 
The impact on the error on 2(3 + 7 by varying the decay model is studied as a function of 
the luminosity and shown in Figure 0] 

From this study the first conclusion is that the present measurements |E],!I2] are n °t 
sufficient to fix the Dalitz model and give a precise indication on the error which can be 
obtained on 2(3 + 7. In fact, the precision on 2(3 + 7 strongly depends on the variation of 
the branching fractions of neutral B into D** states within the measured errors. There 
is also some dependence on the presence or not of non resonant contributions and on the 
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Table 3: Average absolute uncertainty obtained on the parameter 2/3 + 7 as a function of 
the fit configuration for an integrated luminosity of lab' 1 . The configurations are explained 
in the text. See also Figured 
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Figure 5: Average absolute uncertainty obtained on the parameter 2(3 + 7 as a function 
of the fit configuration. The magenta thick stars refers to fit with no background while 
red dots and blue square correspond to fit where the level of background is set at 50% and 
30%. The background has been simulated flat over the Dalitz plot. The fit configuration 
are explained in the text. 



values of strong phases for the most relevant resonances. We also show the impact of the 
variation of r = goZ^"°^o from 0.3 to 0.5. 

The model using only a non resonant contribution for the V u b component is also given 
as a reference since it was taken as example in the feasibility study presented in 
The corresponding error on 2/3 + 7 is a factor 3-6 better with respect the one obtained 
using a set of realistic decay models. This model has not been further considered as it is 
contradiction with present data [Zj. 

3.2 Model Independent fit 

An important feature of this method is the possibility to have a model independent 
determination of 2(3 + 7 by fitting the parameters of the decay model. We consider a 
sample corresponding to 1 ab -1 of collected statistics. The reference error is obtained 
as previously by fitting only 2/3 + 7 an d by fixing all the other parameters. We try to 
progressively release the other parameters which characterise the decay model. 

First we relax the D** and the non-resonant V c b components. We managed to fit 
the three amplitudes and phases and the error on 2(3 + 7 increases by 10%. So far we 
have used only tagged events since they are the only one carrying information on 2(3 + 7. 
Nevertheless also untagged events could play an important role in fitting the decay model 
parameters. In practice we can add the following term to the likelihood given in Eq0] : 

Puntag(B° - D-K°7T + ) = ^k±^k (10) 

The same fit is repeated with untagged events included in the sample. The error on 2/3 + 7 
decreases by about 10%. In the following untagged events are used. 
The following test is to leave free the strong phases of the D** and the non-resonant V u b 
components. In this case the error on 2(3 + 7 increases by about 20%. An extra increase 
of 10% is obtained if the amplitudes of those components are further left free to vary in 
the fit. 

One can also try to estimate the effect of the presence of the background in the Dalitz 
plot. For simplicity we take a flat background over the Dalitz plot, being aware that reality 
could be different. We assume that the discrimination between signal and background 
events is performed using additional variables (e. g. the reconstructed B mass) and we 
fix the signal over background ratio in the CP violation fit. We run our simulation by 
considering a signal over background (S/B) ratio of 30% or 50%. The result is that the 
error on 2(3 + 7 increases by 25% and 50%, respectively. Furthermore if S/B is equal to 
50% it seems difficult to perform a fit where all V c b and V u b D** and the non-resonant 
components are left free in the fit. 

As by-product of this analysis, the Dalitz model can be fitted from the data if the level of 
background is not too large. The V c b D** related amplitudes and phases can be fitted with 
precisions lying between 20-30% . Furthermore both the phases and the amplitudes for 
D** from V u b processes are determined with errors which can go up to 50%. In addition 
the excited K resonnaces (ir o *(1430) ± , ir 2 *(1430) ± , iT*(1680) ± ) which concerns only V cb 
processes can be precisely determined from data with relative uncertainties better than 
25%. 



4 Conclusions and perspectives 



We have performed a full feasibility study for measuring 2(3 + 7 with time dependent 
Dalitz analysis using D~K°tt + final states. For this study a realistic decay model, based 
on the available experimental information has been elaborated. 

We conclude that the error on 2(3 + 7 strongly depends on the decay model and that 
the currently available experimental information is not enough for a precise estimate. By 
varying the different parameters in the currently allowed range, the error on 2(3 + 7 can 
vary within a factor of about 2.5. Only a complete Dalitz analysis could tell us where we 
stand. This is possible on data. 

In fact we have shown that if the level of background is not too large (< 50%) the full 
decay model can be obtained by fitting the Dalitz plot using both tagged and untagged 
events. 

All the amplitudes and strong phases can be fitted and, in case of the Dalitz decay 
model given in Table the precision on 2/3 + 7 could be of about 50% with 500 fb _1 and 
should not be limited by the systematical uncertainties. This implies that the precision 
on 2(3 + 7 could be reduced at few percent level in case this analysis is performed at a 
high luminosity B factory [TT] . 
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